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We demonstrate the possibility of existence of indirect moving excitons in graphene. The binding 
is conditioned by the trigonal warping of conic energy spectrum. The binding energies are found for 
the lowest exciton states. These energies essentially depend on the value and direction of exciton 
momentum and vanish when the exciton momentum tends to the conic points. The ways to observe 
the exciton states are discussed. It is shown that the statistics of graphene at finite temperature 
experiences corrections caused by the exciton states. 

PACS numbers: 71.35.-y,73.22.Lp,73.22.Pr, 78.67.Wj, 65.80.Ck 



Introduction. Graphene is a remarkable material due 
to its unusual conic electron energy spectrum. One of the 
important consequences of this singularity is zero effec- 
tive mass of electrons near the conic point that results in 
a high potential mobility of charge carriers in graphene. 
Just this fact makes graphene very perspective for elec- 
tronics. Another important property of free-suspended 
graphene is potentially large dimensionless interaction 
constant that strongly distorts a single-particle spectra. 

The purpose of the present paper is an envelope- 
approximation study of possibility of the Wannier-Mott 
excitons formation near the conic point in a neutral 
graphene. There were a lot of studies considered the 
excitons in graphene. Unlike the present paper, the term 
" exciton" in all these papers is used for many-body (" ex- 
citonic") effects^, exciton insulator with full spectrum 
reconstruction, or exciton-like singularities originating 
from saddle points (van Hove singularity) of the single- 
particle spectrum^ where the real electron-hole binding 
does not occur. On the contrary, our goal is pair bound 
states of electrons and holes. There is widely accepted 
opinion that zero gap in graphene forbids the Mott exci- 
ton states (see, i.g.^). This statement valid in the conic 
approximation proves to be incorrect beyond this approx- 
imation. It will be shown that exciton existence is condi- 
tioned by the trigonal corrections to the conic spectrum. 

The trigonal warping of the energy spectrum is respon- 
sible for the electron-electron scattering selection rules in 
the neutrality point and, consequently, a specific mecha- 
nism of the valley currents in graphene^. Here we shall 
demonstrate that the presence of the warping drastically 
affects the possibility of electron-hole binding also. 

We consider graphene with weakly interacting charged 
carriers, when the interaction constant g = e 2 /vx <C 1, 
where e is electron charge, v is electron velocity in the 
conic approximation, \ is the half-sum of dielectric con- 
stants of surrounding media, Ti — 1. First, let us dis- 
cuss the behavior of a classical electron-hole pair with 
conic spectra vp. Let electrons and holes have parallel 
momenta. Due to equality of velocities, the distance be- 
tween e and h remains constant. The interaction between 
them changes the momenta but does not change their dis- 
tance. This behavior strongly differs from a usual pair 



with quadratic spectra where the interaction produces a 
well in which e and h oscillate. The quantum counter- 
part of this oscillation is an e-h pairing. Beyond the conic 
approximation the series expansion of the pair kinetic en- 
ergy in the momentum determines the finite inverse mass 
(second derivative of energy with respect to the momen- 
tum) that leads to the variation of the e-h distance. How- 
ever, the oscillations occur only if the sign of this mass is 
positive. Below, we demonstrate that this sign depends 
on the momentum angle, that means the binding exists 
in particular sectors of the momentum space and is for- 
bidden in others. 

Exciton states. We start from the tight-binding elec- 
tron Hamiltonian 
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where 7 = 2u/aV3, 61 = (0, l)a/%/3, b 2 = (-3, -V3)a/6, 
63 = (3, — \/3)a/6, a = 0.246 nm is the lattice constant. 
The spectrum reads as 
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We will consider the case when the binding is weak and 
the exciton states are composed from the narrow group in 
the momentum space. Near points K and K' , where p = 
±K, K = (l,0)47r/3a, the spectrum is conic: e±(p) = 
±v\p — K| and e±(p) = ±f|p + K|, correspondingly. In 
the envelope approximation the two-body Hamiltonian is 

H ex =H e (p e )®I h -I e ®H h (p h )+Ie®IhV(r e -r h ), (2) 

where Hh(p) — —H e (—p) is the hole Hamiltonian, r e ,/i 
and p ej h are electron and hole coordinates and momenta, 
V(r) = — e 2 /vr is the potential of electron-hole interac- 
tion. The electron and hole Hamiltonians relate to the 
electron and hole subspaces of quantum numbers. Hamil- 
tonian ([2]) has a 4-fold set of states, one of which belongs 
to the ordinary exciton, namely, a pair with " an electron 
in the conduction band and a hole in the valence band" . 

Let us carry out the unitary transform U = 
U e (p e )Uh(ph) of the Hamiltonian ([2]), where U e (p e ) and 
U h (ph) diagonalize H e (p e ) and H h (p h ), respectively. 
The smallness of the binding energy means large size 



2 



of exciton wave function and the smallness of the spa- 
tial derivatives; in that case one can neglect the non- 
commutativity of U and V(r e — r^). As a result, we get 
to the Hamiltonian: 



H = e+(p e ) + e+(p h ) + V(r e - r h ). 



(3) 



The electron and hole momenta can be expressed via 
pair q = p e + Ph and relative p = p e — Ph momenta. The 
momenta can be situated near the same (q —¥ k <C 2K) 
or near the opposite conic points (q = 2K + k, k <IC 2K) 
(see Fig. [J). 

Further simplifications are done assuming that the 
binding energy of the pair is less than kinetic energy vk. 
The result is the Schrodinger Hamiltonian 
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where the coordinate system with basis vectors ei = k/fc 
and e2_Lei is chosen. In the conic approximation we have 
m = k/v, M = oo. To obtain finite M, one should go 
beyond the conic approximation. For the indirect ex- 
citon (electron and hole from different cones) with the 
momenta close to z^2K, (v = ±1), we have 1/M = 
uva cos(3^k)/4v / 3, where 0k is an angle between k and 
K. Effective mass M is directly determined by the trig- 
onal contributions to the spectrum and is parametrically 
large: r\ = M/m — A\f?>/{vka cos 30k) ^> f- The sign of 
M is determined by v cos 30k- If v cos 30k > 0, electron 
and hole tend to bind, otherwise to run away from each 
other. Thus, the binding is possible for v cos 30k > 0. 

The similar reasoning for the pair from the same val- 
ley shows that the trigonal contributions to 1/M cancel 
each other. The contributions of a higher order give, 
in this case, a negative mass 1/M: 1/M = — kva 2 (7 — 
cos60k)/32. As a result, the electron- hole binding is for- 
bidden for particles from the same valley. 

The indirect exciton existence needs, positively defi- 
niteness of the effective mass tensor. The corresponding 
condition valid apart from the conic point reads 

(1 + s + t_)(l + s + t+)(l + t- + t+) > 0, (5) 

where s = cosak x , t± = cos((k x ± \JHk y )a/2). 

Fig. Q] shows the domain of exciton existence in the 
momentum space. This domain covers a small part of 
the Brillouin zone. 

The problem of Coulomb states for the Hamiltonian (|4]) 
can be considered using strong anisotropy of the energy 
spectrum. The kinetic energy of the free pair vk is the 
main part of the exciton energy vk — e n N, where £ n jv > 
is the binding energy, quantum numbers n and N numer- 
ate the exciton state. A large ratio of masses 77 3> 1 re- 
minds of the problem of molecular states, where m and M 
represent electron and ion masses, correspondingly. So 
the way of solution is similar to the problem of molecular 
levels. First, we cancel pf/2M, fix the "ion" coordinate 
"x" and determine the energy terms vk — e n (x), then use 
these terms as ions potential energy. The binding energy 
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FIG. 1. Left: Direct and indirect exciton formation in 
terms of single-particle spectrum. Right: Relief of the single- 
electron spectrum. Energy runs from in the conic points to 
37 in the center. The trefoils bound the permitted domains 
of q/2 for indirect exciton momentum q. 



e n N in the N-th "ion" state on the n-th "electron" term 
is determined by the Schrodinger equation for "ions" 



{ P \/2M - e n (x) + e nN ) V nN (x) = 0. 
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The order of magnitude of the lowest levels is the Bohr 
energy e# = me 4 /% 2 . 

Let us replace the potential of electron- hole interaction 



by V(r) 
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\y\) which has the same asymp- 



totics. This approach permits to approximate the energy 
terms at small The replacement converts the 

Schrodinger equation to the analytically solvable one: 
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The solution of Eq. ([7]) satisfying zero boundary condi- 
tion at \y\ — > co (e(x) > 0) is 



i>(y) = c ±e - p ™+H)(|y| + \ x \) 

U[l-(P(x)a B )-\2,2P(x)(\y\ + \x\)] 
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where clb — x/me 2 is the Bohr radius, U[a, b, z] is the 
confluent hypcrgcomctric function, P(x) — i/2me(a;), 
coefficients C± correspond to y > and y < domains. 
Even and odd solutions satisfy the boundary condition 
^eveni®) = and *Podd(0) = 0, accordingly. This gives 
the energy terms at x <C as- eq(x) = 2eslog 2 
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where n > 1. Using the energy terms and Eq. ([5]) we 
obtain the energy levels in the quasiclassical approxima- 
tion: 
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FIG. 2. Left: Binding energy (in eV) of the ground state 
of indirect exciton in graphene versus wave vector in units 
of reciprocal lattice constant. Exciton exists in the sectors 
shown in left figure. Right: Radial sections of right figure at 
fixed angles in degrees (marked). Curves run up to the ends 
of exciton spectrum. 
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The distance between " ions" levels e nj 
less, than the distance between "electron" levels £ n +i,N — 
£n,N by measure of the parameter 1/rj. 

Application of the variational approach to the 
Hamiltonian Q with full variational wave functions 

cxp (— x 2 /a 2 x — y 2 ja 2 y ) and exp (—^Jx 2 /a 2 + y 2 /a 2 ) gives 

£oo = aeslog 2 ?7, where a = 1/n and a = 4/7r 2 , corre- 
spondingly. These values of a do not strongly differ from 
1/2 found in the exactly solvable model (see Eq. (fTTj) ). 

The binding energy of the lowest level £oo essentially 
depends on the angle </>k, while the angular dependence 
of the other levels touches only k-dependent corrections 
to these levels. The corrections have essential-singular 
behaviors as functions of the parameter r\ different for 
odd and even states. 

We have numerically calculated the binding energy of 
the lowest exciton level £oo using the adiabatic approach. 
Figures[2]show the dependencies of £oo on the exciton mo- 
mentum. The calculations have been done by two steps: 
finding the energy terms and solving the Schrbdinger 
equation in the y-direction. 

It should be emphasized that full exciton energy vk — 
£ mn is always positive lying in the electron-hole contin- 
uum. However, the momentum conservation preserves 
the exciton from the decay, unless the scattering pro- 
cesses (weak by assumptions) are taken into account. 

Statistics of graphene with excitons. Another unusual 
property of graphene is participation of excitons in equi- 
librium thermodynamics. The presence of excitonic 
states changes statistics and thermodynamic properties 
of graphene at the finite temperature T . Excitons are 
Bose-particles and the number of equilibrium excitons 
is indeterminate. Hence, in the Dirac point, they have 
a zero chemical potential. Excitons in graphene do not 
Bose-condensate due to their gaplcss spectrum. The exci- 
ton contribution to extensive thermodynamic parameters 
can be calculated by subtracting the contribution of free 
pair (that is included into the electron and hole contribu- 



tion) from the bound pair. The correction to the energy 
is 



SE 



r 2 



d 2 p 



vp - e nN 



up 



1 



1 



(14) 



where k B is the Boltzmann constant. Taking into account 
the smallness of e n N <C vp, we find the correction to the 
specific heat 

C ex « a{gk B T/v) 2 [ Cl log 2 6/T - c 2 log 6/T + c 3 ],(15) 

where 9 = A^v/k B a, c x « 9.18, c 2 ~ 3.46, c 3 w 8.93. 

The contribution of bound pairs should be com- 
pared to that of the free electrons and holes C e +h — 
3 Cl (k B T/2v) 2 . 

The specific heat of graphene at low temperature is 
mostly determined by the acoustic and softer vertical 
phonon modes. The acoustic phonons give contribution 
C p h oc T 2 /v 2 , where v s is the sound velocity. This should 
be compared to the electrons and holes contribution in 
the Dirac point cx T 2 /v 2 . 

How to observe excitons in graphene? Let us discuss 
the possibility of observation of the indirect excitons in 
graphene. In a semiconductor with the energy gap ex- 
citons contribute to the absorption and, especially, the 
emission of light. The absence of the energy gap makes 
observation of excitons in graphene a different problem 
from the case of usual semiconductor, because exciton 
energies are distributed between zero and several tenth 
of eV that smears up the exciton resonance. The large 
momentum of indirect exciton blocks both direct optical 
excitation and recombination. However, the slow recom- 
bination and inter- valley relaxation preserve the excitons 
(when generated someway) from the recombination or 
the decay. The frequency of photons with a small wave 
vectors, which can produce excitons, is very low. In that 
case, the solution can be found using a structure where 
electrons obtain an additional in-plane momentum from 
phonons with wave vectors K, impurities, the artificial 
grating or the naturally riffle of the graphene plane. 

Phonon-assisted indirect exciton optical absorption 
has thresholds determined by one of phonon frequencies 
f2,-(K), where i numerates the phonon type. Correspond- 
ing contribution to the absorption near the threshold 
fii(K) is determined by the exciton density of states and 
is proportional to (u> — Qi)9(uj — fl;), where uj is light 
frequency. 

The other way perspective to observe excitons is based 
on a monolayer character of graphene. Being placed in 
vacuum, the graphene layer can serve as electron tube 
grid. This permits to study the scattering of vacuum 
electrons (forward or backward) on the graphene. The 
dependence of the energy loss maximum on the scatter- 
ing angle gives the spectrum of the Bose excitation, in 
particular, excitons. To provide the necessary momen- 
tum transfer, the external electrons should have energy 
~ lOeV. The energy resolution of modern methods of 
low-energy electrons scattering is sufficient to resolve ex- 
citons with energies of some meV. 
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The electron subsystem specific heat can be isolated 
in dynamic experiments with quick electron heating by 
means of pulse voltage or light. In such conditions the 
electron temperature is separated from the phonon one 
and the electrical or optical responses will reflect the time 
evolution of the electron temperature connected with 
the electron cooling. Assuming quick establishment of 
the temperature equilibrium between electrons and exci- 
tons and slow establishment of the equilibrium between 
them and phonons one can find the the ratio of tem- 
perature changes just after pulse 5T and after equaliz- 
ing of phonon and electron subsystem temperatures ST^ : 
STo/ST^ = 1 + C ph /(C e+h + C ex ). This gives the ther- 
modynamic way to observe excitons in graphene. 

Discussion and conclusions. In the above considera- 
tion the many-body corrections have been neglected. Ac- 
tually, the e-e interaction leads to the electron self-energy 
gvk\n(l/ka), g = g/A- # 2 (5/6 - In 2) + 0(g 3 ) 6 . We have 
neglected these corrections, assuming that not only g, 
but typical gh\(l/pa) is less than unity. The presence 
of this logarithmic term strongly changes the behavior of 
the single-electron spectrum at small k and gives a neg- 
ative contribution to M. However, this does not forbid 
the exciton absolutely. In fact, the interaction-induced 
correction to the inverse mass l/m.;„ t — ~gv/k should 
be added to the warping correction. In a particular case 
of small interaction constant g, mass M becomes nega- 
tive at low k and near the boundaries of the existence 



sectors. This narrows the range of exciton existence and 
shifts it from the conic point. 

Another remark concerns the renormalization of the 
dielectric constant x, caused by the vacuum polarization 
of graphene!. The renormalization reduces the Coulomb 
interaction at large distance, as if the dielectric constant 
of external medium becomes logarithmically large. This 
fact improves the applicability of the weak-interaction 
model, even in the case of free-suspended graphene. 

The excitons considered in the present paper arise near 
the conic points, unlike the saddle-point exciton states 
with energies about 5 eV studied io^r^r—. The latter, in 
our opinion, correspond to the e-h scattering resonance, 
rather than bound states. The other investigations of ex- 
citon resonances is based on the electron spectrum recon- 
struction caused by large interaction constan t 14 ' 15 . How- 
ever, the possibility of gapless excitons in a system with 
small g considered here had not been studied earlier. 

In conclusion, we have demonstrated the existence 
of moving indirect excitons in monolayer graphene and 
found their energies. The excitons are conditioned by 
the presence of the warping of the electron spectrum. 
The excitons spectrum has no gap. The many-body in- 
teraction corrections to the single-particle spectrum, in 
concurrence with warping, lead to the destruction of ex- 
citons with small wave vectors. However, the excitons 
survive if the interaction constant is small. 
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